The dynamics of vortex lines in the three-dimensional complex Ginzburg-Landau 
equation: instability, stretching, entanglement, and helices 
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The dynamics of curved vortex filaments is studied analytically and numerically in the framework 
of a three-dimensional complex Ginzburg-Landau equation (CGLE). It is shown that a straight 
vortex line is unstable with respect to spontaneous stretching and bending in a substantial range of 
parameters of the CGLE, resulting in formation of persistent entangled vortex configurations. The 
boundary of the three-dimensional instability in parameter space is determined. Near the stability 
boundary, the supercritical saturation of the instability is found, resulting in the formation of stable 
helicoidal vortices. 
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I. INTRODUCTION 

Analysis of " universal" models plays a central role in 
contemporary nonlinear dynamics. Such models, as with 
Ginzburg-Landau and Swift-Hohenberg equations, allow 
for quantitative description of arbitrary nonlinear sys- 
tem at the threshold of specific instabilities. The com- 
plex Ginzburg-Landau equation (CGLE), derived some 
20 years ago by Newell and Kuramoto ||] has be- 
come a paradigm model for a qualitative description of 
weakly-nonlinear oscillatory media (see for review [^). 
Under appropriate scaling of the physical variables, the 
equation assumes the universal form 



dtA = A-{l + ic)\AfA + (1 + tb)AA, 



(1) 



where A is the complex amplitude, b and c are real pa- 
rameters, and A = + dy + is a three-dimensional 
Laplace operator. The parameter b is the ratio of dis- 
persion to diffusion and c is the ratio of conservative to 
dissipativc nonlinearity. Although the equation is for- 
mally valid only at the threshold of a supercritical Hopf 
bifurcation, it has been found that the CGLE often repro- 
duces qualitatively correct phenomenology over a much 
wider range of the parameters. As a results, the predic- 
tions drawn from the analysis of the CGLE (mostly in 
one and two spatial dimensions, see e.g. ^-0]) were re- 
cently successfully confirmed by experiments in optical 
and chemical systems [||J^]. Moreover, some results ob- 
tained from the CGLE (for example, symmetry breaking 
of spiral pairs) was instructive for interpretation of ex- 
periments in far more complicated systems of chemical 
waves ijl^ and colonies of amoebae ]ll| , p^ 

In three dimensions the point singularity at the center 
of the spiral becomes a line singularity, known as a scroll, 
or vortex filament. The filaments can be open (scrolls), 
closed (vortex loops and rings), knotted or even inter- 
linked or entangled. Scroll vortices had been observed 



in slime mold (l3| , heart tissues [|l4|, gel- immobilized 
Belousov-Zhabotinskii (BZ) reaction [[15[. Long-lived en- 
tangled vortex patterns in three-dimensional BZ reac- 
tions were observed recently by the group of Winfree us- 
ing advanced optical tomography technique Com- 
plicated vortex configurations have also been observed 
in numerical simulations of reaction-dilfusion equations 

Theoretical investigation of scroll vortices in reaction- 
diffusion systems was initiated by Keener j20j, who de- 
rived the equation of motion for the filament axis. In 
particular, he obtained that the collapse rate is propor- 
tional to the local curvature of the filament, leading to a 
collapse of a vortex ring in a finite time. The existence 
of non- vanishing vortex configurations and expansion of 
vortex loops, observed also in numerical simulations of 
reaction-diffusion equations, was associated with "nega- 
tive line tension" of the vortex filament . 

Recently, the dynamics of three-dimensional (3D) vor- 
tex lines in the CGLE has attracted substantial attention 
pl|-p3| . As a definition of a vortex line we accept a line 
singularity of the phase of a complex function A. Gabbay 
Ott and Guzdar applied a generalisation of Keener's 
method for a scroll vortex in reaction-diffusion systems 
po| . They derived that the ring of a radius R collapses 
in finite time according to the following evolution law 



dR 



R 



(2) 



In addition, there is no (at least, in the first order in 1/ R) 
overall drift of the vortex ring in the direction perpen- 
dicular to the collapse motion. The collapse rate (called 
often " line tension" ) v = 1 -\-b^ appears to be in a rea- 
sonable agreement with the simulations |^^. This result 
generalises Keener's ansatz by including the curvature- 
induced shift of the filament's wavenumber. Thereby, as 
follows from Eq. (0), vortex loops initially existing in the 
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system will always shrink (if, of course, there is no bulk 
instability of the waves emitted by the vortex filament), 
and under no condition can the vortex loop expand. 

In this article we show that under very general condi- 
tions, and in an extensive part of the parameter space 
vortex filaments expand and bend spontaneously and re- 
sult in persistent vortex configurations even if there is 
no bulk instability of the emitted waves and the spiral 
wave is stable in a two-dimensional system. A prelim- 
inary account of this work had been published in Ref. 
p3[ . Wc have shown that vortex loops may expand for 
any value of b above some critical value bc{c). The criti- 
cal value 6c (c) can be relatively small for not too large c. 
For example, our analysis predicts that bnK, 2 for c — Q. 
We prove that for b > be the equation (g) is not valid, 
because formally higher-order, but in fact singular cor- 
rections, omitted in Eq. cause severe instability of 
the filament and persistent stretching. This instability is 
a three-dimensional manifestation of the two-dimensional 
core instability of spiral waves (called in Ref. |^] acceler- 
ation instability). Its origin can be traced back to the 
breakdown of the Galilean invariance of the CGLE for 
any e = 1/b ^ 0, causing spontaneous acceleration of the 
spiral waves Whereas in 2D the instability is rela- 
tively weak, the situation is different in three dimensions. 
Using combined analytical and numerical methods, we 
have proven that three-dimensional instability of the vor- 
tex filaments persists far beyond the core instability limit 
of two-dimensional spiral wave and typically has a much 
higher growth rate. It ceases to exist only when the core 
modes becomes strongly damped. This instability is not 
driven by "negative line tension'. It develops from a non- 
trivial response of the filament core to bending, which re- 
sults in additional " acceleration" terms in the equation of 
filament motion. As we will show, the bending of the fil- 
ament greatly enhances the instability, and may result in 
formation, after some transient, of a highly entangled and 
dense vortex configuration. The "high dispersion limit" 
6 S> 1 is readily fuUfilled for many physical and chem- 
ical systems (actually b larger than about 2 is enough 
to have the effect). For example, in the context of non- 
linear optics, where the CGLE can be derived from the 
Maxwell-Bloch equation in the good cavity limit |Q , this 
parameter is very small: e = 1/b ^ 10^^ — 10^^. For an 
oscillating chemical reaction the diffusion rates of various 
reacting components can be varied over a wide range by 
adding extra chemicals. 

The structure of the article is following. In Sec. II 
we present an analysis of the filament motion in the 
"high dispersion" limit e = 1/6 <C 1. This allows us 
to prove instability of a straight vortex filament with 
respect to bending and stretching on the basis of a 
computer-assisted analytical procedure. In Sec. Ill we 
present the results of numerical simulations of the three- 
dimensional CGLE. We discuss the properties of a spatio- 
temporal intermittency which we have found in our sim- 
ulations. In Sec. IV we consider the formal stability 
analysis of a straight vortex filament with respect to pe- 



riodic perturbations along the filament axes. Using a nu- 
merical matching-shooting algorithm, we have calculated 
the spectrum of eigenvalues and determined the stability 
limit of a three-dimensional vortex line in the 6, c plane. 
In Sec. V we present a weakly-nonlinear analysis for the 
vortex filament near the three-dimensional stability limit 
showing the existence of traveling helix solutions. After 
some concuding remarks in Sec. VI we present in the Ap- 
pendix an analytical derivation of the equation of motion 
of the vortex line and the stability boundary in the limit 
of the perturbed nonlinear Schrodinger equation (NSE). 

II. INSTABILITY OF WEAKLY-CURVED 
FILAMENTS IN A HIGH DISPERSION LIMIT 

A. Derivation of the equation of motion 

As a test for instability, we consider the dynamics of 
a weakly curved vortex filament in the high-dispersion 
limit 6 3> 1. We apply the generalisation of the method 
of Ref. 1^ for the case of 3D vortices, and make pertur- 
bations near the 2D spiral wave solution to the CGLE. 
For convenience, we redefine r r/Vb. Then Eq. (|l|) 
assumes the form 

dtA^ A- {l + ic)\A\^A + {e + i)AA. (3) 

In the following discussion in this section, we assume < 
e <^ 1 to be a small parameter. Our objective is to 
relate the acceleration of the vortex filament dtv with 
the velocity v and local curvature k of the filament. 

In order to develop the perturbation theory for a 
weakly-curved vortex filament in 3D, we begin with the 
stationary one-armed isolated spiral solution to Eqs. 
(|l|;^, which is of the form 

Ao{r, 9) = F{r) expi[ujt ±9 + V'(r)] , (4) 

where (r, 9) are polar coordinates, lu — —c — fco(l ~ ^c) is 
the rotation frequency, and fco is an asymptotic wavenum- 
ber. The real functions F and tp have the following 
asymptotic behavior F{r) \J'v^~(}i^ ,V''('') ^ 
for r ^ oo and F(r) ~ r yip'i''') ~ ^ for r 0. The 
wavenumber fco of the waves emitted by the spiral is de- 
termined uniquely for given e, c [ p5| . For e = one has a 
type of Galilean invariance and then, in addition to the 
stationary spiral, there exists a family of spirals moving 
with arbitrary constant velocity v — {vx,Vy) 

A{r,t) ^ F{r')e:ipi[uj't + 9 + i;{r') + ^^], (5) 

where r' = r — vt, ut' — u; ~ and the functions 

F,4> are those of Eq. (^). (This invariance holds for 
any stationary solution). For e 7^ the diffusion term 
~ eAA destroys the family and leads to slow acceler- 
ation/deceleration of the spiral proportional to eu, de- 
pending on the value of e. As was found in Ref. the 
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equations of motion of the spiral core for e ^ 1 assume 
the form 



dtv + eKv = 
with the "friction" tensor 

7>- ( ^xx ^xy 

K K 



(6) 



(7) 



dtv + exv = 0, 



(8) 



where v = Vy + iv^ is a "complex" velocity and x — 
Kxx + iKxy is a "complex" friction coefficient. It was 
shown in the Ref. that the coefficient K^x < for 
e ^ 1, which implies instability of the spiral core (see 
Fig. |l|). In the hmit of the weakly-perturbed NSE (c 
oo) the friction x is of the form (see Appendix A): x ~ 
— 3.45c'^ exp[— ctt] (1 ± i3.92). Since in general K^y ^ 0, 
the spiral core moves on a logarithmic spiral trajectory. 




small. We assume that the velocity vector v lies in the 
plane perpendicular to the vector tangential to the vortex 
line and may in general depend (slowly) on the arclength 
s. Substituting ly into Eq. (|), we obtain at first order in 
e the linear inhomogeneous (assuming dtv,v, k <^ 1 and 
neglecting higher-order corrections dsV,d^v etc.) equa- 
tion: 



Because of isotropy the elements of the tensor satisfy the 
condition: K^x = Kyy.K^y = -Kyx- Eq. (||) can be where 
written in a more compact form: 



L{W,W*) = H, 



(10) 



L{W, W* 



iAW 

^{l + ic)F'^{W + W*) 
j_ 1 a2 



(11) 



and A ^ + ^dr + -^d^ - -^dr{rdrF). The inhomo- 
geneity H is of the form 

H = ^dt-vF - idtipoF - ev {{VF + iFV{e + ?/>)) 
+ K{e + i){diF + iFdi{e + 4,)), (12) 

where (r, 9) are polar coordinates in a local plane spanned 
on the vectors of normal and binormal, and we used the 
expansion for the Laplace operator in the local basis in 
the limit of small curvature k and torsion t: A — —Kdi + 

+ d~ + dg + .... We use the notation v — vb + ivn, 
where vb , vn are the bi-normal and normal components 
of the velocity, respectively (note that the coordinate x 
is directed along the vector of normal). 

Separating real and imaginary part of W, and repre- 
senting it in the form of a Fourier series 



I'ReW 
\lmW 



E 



An{r) 
Bn{r) 



exp(m0), 



(13) 



FIG. 1. Kxx (solid line) and K^y (dashed line) as functions 
of c. Lines are theoretical results, symbols are the results of 
two-dimensional simulations from Ref. for e = 0.025. 

Let us now consider the dynamics of an almost straight 
vortex line. The analysis is conveniently performed in a 
filament-based coordinate system (see for details [p2|). 
The position in space X is represented by local coordi- 
nates s,x,y^ where s is the arclength along the filament, 
and X = R(s) + iN(s) -I- yB(s), where R is the coor- 
dinate of the filament, N(s) the normal and B(s) the 
binormal vector (see also Fig. 1 in Ref. |2^). In this 
basis the weakly curved filament moving with velocity v 
can be written in the form 



A{r,t) = {F{r') + W{r',9,s)) 



exp 



i (iu't + e + ^0 + + 



r • V 



(9) 



Here r' = r — vt, 0o is a constant phase, and W is the 
perturbation to the spiral solution which we require to be 



we arrive at the set of ordinary differential equations for 
each azimuthal modes yl„,_B„. Since instability occurs 
only for the first azimuthal mode, we consider the equa- 
tion only for n — I (the equation for n — —1 is obtained 
by complex conjugation). At first order in e the inhomo- 
geneous linear equation for the corrections Ai , Bi is of 
the form 



A^i - 2{cF^Ai + i>'F 



d El ^ iBi , 
dr F 



ev 



F' 



rdtv 



F 



2i 4i 
ABi + 2{F^Ai +^'F 



2^' + y(^^' 

d Ax iAi 
dr F r 



2 ) 



iF^ 



2i 



(i^V' + — ) 



4/' 



'-{F^/y 



r ' 



(14) 



where A — 



- (d^F + l/rdrF)/F and 



primes denote differentiation with respect to r. Equation 
( |T^ ) can be formally simplified by the following transfor- 
mation Bi = Bi - erF/4:,Ai = Ai p^]. After simple 
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algebra, Eqs. (|TJ) assume the form (we omit tildes on 

or t r^ 

2i Ai 

or t r^ 

iF ^ 



where r = (Xq^Yq). Using 9tr = v, from Eq. ( |17| ) we 
then obtain: 



(15) 



where v = v — in/e. The solvability condition for Eqs. 
( p^ implies an unique relation between dtv and v. If 
the solvability condition is fulfilled, the solution Ai,Bi 
remains finite at r = and does not grow exponentially 
for r — !■ cx). Slow, power-like growth of the solutions is 
permitted, since ther.h.s. of Eqs (|l|) grows linearly with 
r. Thus the transformation Bi ^ Bi — erF/A does not 
change the solvability condition. Remarkably, Eqs. (|l5| ) 
coincide with the equations analysed in Ref |^ for the 
case of the acceleration instability of a spiral wave in two 
dimensions. The only difference is that the normal veloc- 
ity is modified by the curvature k. Obviously, it results 
in the same solvability conditions as the corresponding 
equations of the two-dimensional case: 



dtV -\- exv = 



(16) 



with the same "complex" friction coefficient x- The equa- 
tion of motion 1l6f) can be written in the matrix form 



(17) 



(9tv + k[ev - kN] = 0. 



Note, that dropping the acceleration term in the equation 
( |l7| ) we recover the result of Ref. |2|] for 6 — + oo, since 
for the ring vn = dtR, k = —1/R. Restoring the original 
scaling r r/Vb, we obtain dtR — ~b^/R. However, 
since in three dimensions the local velocity in general 
varies along the vortex line, even small acceleration may 
cause severe instability of the vortex line, because the lo- 
cal curvature becomes very large. Moreover, deviation of 
the local velocity from the direction of normal will lead 
to stretching and bending of the vortex line. Thus the 
acceleration term, which formally can be considered as 
a higher order correction to the equation of the motion, 
may play a pivoting role in the dynamics of a vortex 
filament. Our subsequent analysis and numerical simu- 
lations verify this hypothesis. 



B. Almost straight vortex 

Let us consider an almost straight vortex parallel to the 
axis z. We can parameterise the position along the vortex 
line by the z coordinate: {X(){z),Y(){z)). Since in this 
limit the arclength s is close to z, the curvature correction 
to the velocity kN is simply kN = {d^Xo,d^Yo) = 9^r, 



dtv + k[ev - a^r] = 0, 



(18) 



Let us now consider perturbations of the vortex that are 
periodic along z. Due to the linearity of the problem, the 
solution can be written in the form r ~ exp[zA:z -I- X{k)t\, 
where A is the growth rate. We immediately obtain the 
following relation for A: 



A' + x(eA-}-fc2) =0. 



(19) 



Let us consider separately two cases: <C e and k ^ e. 
For fc <C e from Eq. (|is| ) we derive that A = — ex + 0(fc^). 
For /c 3> e we obtain A ~ zb^/— (^xTiiKj^fc. There 
always exists a root with a large positive real part: 
A ~ fc ^ e. Therefore, for finite k, the growth rate 
A(/c) may significantly exceed the increment of the ac- 
celeration instability in 2D (corresponding to fc = 0): 
A = —ex = —e{Kxy±iK^y). We can expect that, as a re- 
sult of such an instability, highly-curved vortex filaments 
will be formed. Hence, the "small-curvature" approxima- 
tion considered above can be valid only for finite time. 
Moreover, one may not expect this instability to saturate 
always in a steady-state configuration with finite curva- 
ture (although we will show this possibility to exist). In 
contrast, we suggest that frequently reconnection of var- 
ious parts of the filaments, formation of vortex loops etc, 
will result in persistent spatio-temporal dynamics of a 
highly-entangled vortex state. We expect a fall off of the 
growth rate A at large k which is not captured by the 
small- curvature approximation used her. This effect will 
be included in the treatment below. 



III. NUMERICAL RESULTS 

In order to verify our results numerically and to follow 
later development of the instability, we performed sim- 
ulations of the three-dimensional CGLE. We studied a 
system of 50'^ — 60^ dimensionless units of Eq. (|^) with 
either no-flux or periodic boundary conditions. The nu- 
merical solution was implemented on a parallel 16 pro- 
cessors Origin 2000 computer of the High Performance 
Computing Center at Argonne National Laboratory. We 
applied an implicit Crank-Nicholson algorithm based on 
iteration scheme for inversion of a band matrix. The 
number of grid points was 100'^ — 128'^. We performed 
simulations in the parameter regime away from ampli- 
tude turbulence in two dimensions |6j for various values 
of e and c. 

We numerically verified our theoretical result for the 
growth-rate of linear perturbation for a straight vortex 
Eq. (p^. As an initial condition we selected a straight 
vortex line with small periodic modulation along the z- 
axis. We have measured numerically the growth-rate as 
a function of modulation wavenumber k. The results of 
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our simulations are summarized in Fig. As we see 
from the figure, the growth rate indeed increases initially 
with k, and then falls off for large k. The theoretical ex- 
pression ( |l9|) shows reasonable agreement with the simu- 
lations results for small enough k. The maximum growth 
rate reached for intermediate values of k exceeds in this 
case the growth-rate of the two-dimensional acceleration 
instability (fc = 0) by more then two orders of magnitude. 




FIG. 2. The growth rate ReX{k) as function of k for 
e = 0.02, c = .1. Solid line is the theoretical result for fc <C 1, 
dashed line with symbols the result of numerical solution of 
3D CGLE. Inset: blow-up of small k region. 



a 






c 


d 





FIG. 3. Instability of a straight vortex filament. 3D iso- 
surfaces of \A{x,y, z)| = 0.1 for e — 0.02, c — —0.03, shown at 
four times: 50 (a), 150 (b), 250 (c), 500(d). Similar dynamics 
is observed also for larger value of e. 



6000 




FIG. 4. The dependence of filament length L on time. 
Solid line corresponds to e = 0.02, c — —0.03; dashed line 
corresponds to e = 0.02, c = —0.5. 

The long-time evolution of a straight vortex is shown 
in Fig. ^. As initial condition we used a straight vor- 
tex perturbed by a small broad-band noise. As we can 
see from the figure, the length of the vortex line grows. 
The dynamics seems to be very rapidly varying in time, 
and the line intersects itself many times forming nu- 
merous vortex loops. The long-time dynamics shows, 
however, a saturation when a highly-entangled vortex 
state is developed and the total length of the line can- 
not grow further due to a repulsive interaction between 
closely packed line segments. The dependence of the line 
length on time is shown in Fig. y. As a measure of 
the filament length L we used the following quantity: 
L w Sq^ J Q{Ao — \A{x,y, z)\)dxdydz, where Q is the 
step function: 0(a;) = 1 if x > and = otherwise. 
Aq = 0.1 was used as a threshold value to identify the 
vortex. 5*0 is a constant determined from the condition 
that for the straight line the above integral coincides with 
the actual length. We can identify two distinct stages of 
the dynamics: first, fast growth of the length; second, 
oscillations of the line's length around some mean value. 




FIG. 5. Two snapshots of 3D isosmfaccs of jj4| 

taken in the regime of spatio-temporal intermittency, 
e = 0.02, c — —0.5. (a) Left image corresponds to t « 620 
for Fig. H (b) right image corresponds to t ~ 740. 
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Strikingly, for small enough e, we have observed two 
distinct behaviors of the total vortex length depending 
on the value of c. Above a critical value Cc correspond- 
ing approximately to the so-called convective instability 
range of the two-dimensional spiral (cc — > for e —^ 0), 
the total length approaches some equilibrium value and 
does not exhibit significant fluctuations. On the contrary, 
for c < Cc, the total length exhibits large non-decaying 
intermittent fluctuations around the mean value. Figure 
^ demonstrates the two cases. In Fig. [s] we present snap- 
shots illustrating the structure of the vortex field in the 
intermittent case corresponding to the maximum of the 
length (Fig. ||a) and the minimum (Fig. ||b), respec- 
tively. One sees, that in this situation some segments of 
vortex lines start to expand spontaneously, pushing away 
other vortex fllaments and in such a way making sub- 
stantial vortex free holes around them. Then the insta- 
bility takes over and destroys these almost-straight seg- 
ments of filament, bringing the system back to a highly- 
chaotic state. This dynamics can be considered as a 
three-dimensional spatio-temporal vortex intermittency, 
which is an extension of spiral intermittency discovered in 
the context of the two-dimensional CGLE in Ref. For 
even smaller values of c < — 1 we observed the transition 
to a highly chaotic state, which is an analog of "defect 
turbulence" in the two-dimensional CGLE. In this regime 
small vortex loops nucleate and annihilate spontaneously, 
and large vortex fllaments play no role in the dynamics. 



the collapse motion. However, we have also found that 
small enough rings typically collapse, since then the in- 
stability described above does not have time to develop 
substantial distortions of the ring. Even in this situation 
the ring exhibits a few oscillations of the radius. Proba- 
bly, the oscillations in the collapse rate in Fig. 2 of Ref. 
p2[ are caused by the effect of the acceleration term in 
the equation of motion, which plays an important role 
even in the stable case. 



IV. LIMITS OF THREE-DIMENSIONAL 
INSTABILITY 

The previous analysis indicates instability of vortex 
lines in the limit e — *■ for all c. However, it cannot 
describe the boundary of the instability in the c, e plane 
since the transition to a stable regime occurs for some 
finite value of e. In order to obtain the stability limit 
we performed a full linear stability analysis of a straight 
vortex solution, which is not limited to small k and e. 

The perturbative solution is of the form: 

A = (F(r) + W{r, 6, z, t)) exp [i {cut + 6 + V'(r))] . (20) 



Substituting the ansatz (|20|) into the CGLE and perform- 
ing the linearisation with respect to W as before, then 
separating the real and imaginary part of W and repre- 
senting the solution in the form 



t=100 




t=120 



t=140 




FIG. 6. Sequence of snapshots demonstrating the evolu- 
tion of a vortex ring for e = 0.2 and c — 0.2. 

The evolution of a closed vortex loop is shown in Fig. 
^. Our simulations show that the three-dimensional in- 
stability may prevent the ring from collapse, causing the 
stretching of the loop in the direction transversal to the 



ReW 
IraW 
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E 



Bn{r) 



exp {inO + ikz + X{k)t) , (21) 



we obtain an eigenvalue problem for the eigenvalue \{k). 
Again, we restrict ourself by the analysis of most danger- 
ous perturbation harmonics with n = 1. The resulting 
equations are of the from (compare Eqs. (p^.pi 



KAi - k^Ai - 2(ciF^Ai -I- %l:'F 



d Bi iBi, 



m 

l + e2 



(eAi+Bi) 



or F 



m 

l + e2 



(eSi-Ai), 



(22) 



where ci = (e + c)/ + and C2 = {1 - ec) / {1 + e"^) . The 
functions Ai, Bi are subject to the boundary conditions: 
Ai , Bi are bounded at r = and decay exponentially for 
r — > oo. 

We solved Eqs. ( p^ ) in the range of wavenumbers k nu- 
merically using a matching-shooting method with New- 
ton iterations from NAG, routine d02agf. Since Eqs. ( |2^ ) 
are singular at r = and the solution is required on an 
infinite interval, we applied the following method of so- 
lution of this rather difficult eigenvalue problem. The 
functions Ai,Bi were expanded in a series for r — > 
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and we used the asymptotic expansion for r — > oo. We 
replaced the boundary conditions at r = and r — > oo 
by new boundary conditions at sufficiently small tq and 
sufficiently large r,,. The boundary values were obtained 
from the corresponding asymptotic expansions, while the 
unknown parameters of these expansions were included 
in the shooting-matching procedure. Thus, the numeri- 
cal matching procedure was applied on a finite interval 
To < r < Tg, where tq was typically 10~^ and was 
gradually increased until the eigenvalue A approaches its 
asymptotic value. Since the unperturbed functions F, tp 
are known only numerically, we determined them in the 
same matching routine, solving three nonlinear equations 
(for F, F' and tp' ) and eight first-order linear equations 
(the functions Ai,Bi are complex). We typically used 
5000 mesh points on the interval of integration. 

The spectrum of ReX(k) for two values of e is shown 
in Fig. 0. A typical localized core eigenmode, corre- 
sponding to this three-dimensional instability is shown 
in Fig. H As we can see from Fig. |^, ReX{k) indeed falls 
off for large k. As expected from the previous analy- 
sis, the three-dimensional instability persists beyond the 
two-dimensional core instability. From Fig. ^ we see that 
for e smaller than about 0.3 one has i?eA(0) < (core is 
stable in two dimensions), however there is an instability 
for finite k. The relation between the eigenvalue problem 
Eq. ( p^ ) and the acceleration instability is presented in 
Appendix B. 

We systematically tracked the boundary for the three- 
dimensional instability from the condition maxi?eA = 0. 
The results are shown in Fig. ^. As one can see, the 
three-dimensional instability occurs in a much wider pa- 
rameter range then the two-dimensional core instability. 
Moreover, the typical growth-rate in three dimensions is 
much higher than in two-dimensions. 

We expect that for c — oo the critical line of the 
three-dimensional instability approaches again e = line, 
similarly to the two-dimensional core instability line 
However, our numerical matching procedure failed to 
converge for this region of the parameters since the crit- 
ical mode becomes less localized for c ^ — oo. 
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FIG. 7. ReX and ImX as functions of k for e = 0.3 (solid 
line) and e — 0.14 (dashed line) for c — 0.5. 
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FIG. 8. ReAi,ImAi (solid lines) and ReBi,ImBi 
(dashed lines) as functions of r, obtained from numerical so- 
lution of Eqs.( H) for e = 0.3, c = 0.5 and k = 0.2, corre- 
sponding to A = 0.0163 + i0.122. 

V. WEAKLY-NONLINEAR ANALYSIS 

Our numerical simulations near the stability bound- 
ary of the three-dimensional instability have revealed a 
striking result: the instability saturates, leading to the 
formation of a stable travelling helix solution (see Fig. 
[To| ). The pitch of the helix is determined by the most 
unstable wavelength, and the radius of the helix vanishes 
when approaching the stability boundary. The existence 
of the stable helix solution can be considered a result of 
saturation of the Hopf bifurcation which is thus of super- 
critical nature. 




e 



FIG. 9. Stability limits in three-dimensions (solid 
line) and two dimensions (dashed line), obtained from lin- 
ear stability analysis. Symbols represents the limit of the 
two-dimensional instability, obtained by direct numerical sim- 
ulation of the OGLE from the Ref. Q. Vortex lines and 
two-dimensional spirals are stable to the left of the respec- 
tive lines. 
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FIG. 10. (a) Stable travelling helix solution, obtained 
numerically for e = 0.3 and c — 0.5. (b) The coordinates 
of the vortex line x{t) (solid line) an y{t) (dashed line) as 
functions of time for some fixed z. A supercritical character 
of the bifurcation is apparent. 

Although we have never observed saturation of the ac- 
celeration instability in two dimensions, here the situ- 
ation is different, since the most unstable mode has a 
different spatial structure then the acceleration mode in 
two dimensions. The helix near the stability limit can be 
described in the framework of a weakly-nonlinear analy- 
sis for the relevant order parameter, which characterises 
the local helicity of the vortex line. 

The structure of the amplitude equations can be un- 
covered from the following consideration. The growth- 
rate A(fc) is symmetric with respect to k: A(fc) = A(— fc). 
Therefore, at the transition point two modes with kc and 
— fee, corresponding to counter-propagating waves of op- 
posite helicity, simultaneously become unstable. In gen- 
eral, n{k) — Im\(k) ^ 0. Thus, at lowest order we shall 
obtain two coupled one-dimensional Ginzburg-Landau 
equations for counter-propagating waves with the am- 
pHtudes U and V . Close to the threshold the structure 
of these equations becomes universal: 

dtU + vpd.u = \{k,)u + ]^dl\[k,)dlu 
dtV - vpd^v = \(k,)v + ]^dlx{k,)dlv 



{ai\V\^ + a2\U\^)V 



(23) 



where Vp = dk^{kc) and ai,a2 are complex constants. 
The coefficients of the linear problem can be determined 
from the solution of the eigenvalue problem Eqs. (22). 
To determine coupling coefficients ai_2 some additional 
analysis is required. In principle, they can be extracted 
from three-dimensional simulations. 

The system ( |2^ ) has been studied intensively in vari- 
ous physics contexts (see e.g. [|8-30 ). It exhibits diverse 
types of dynamic behaviors, including travelling waves, 
domain walls, sinks, shocks, and spatio-temporal chaos. 
Thus, in our interpretation, this chaotic behavior corre- 



sponds to slow nonperiodic spatio-temporal deformations 
of the helicoidal structure of the vortex line. 

Equations (|2^) are simplified drastically if one of the 
counter-propagating waves is suppressed, which is the 
case for Reai > Rea2- Then this system is reduced to a 
single one-dimensional complex Ginzburg-Landau equa- 
tion, which is of the form (in a moving frame) : 

dtU^X{kc)U+^dl\d'^,U~ai\U\^U (24) 

We estimated the parameters of Eq. (|24| ) from our linear 
analysis and simulations (see Fig. |1C| b). We obtained 
A(/cc) ~ 0.02095 -I- i0.1431, dl\{k^)/2 = 0.311 - i0.5315 
and fli — ag(l — iO.3675), where is a parameter 
which can be scaled out. For this set of parameters 
of the one-dimensional CGLE (c — —0.3675 and h = 
—0.5315/0.311 w —1.708) the homogeneous solution to 
Eq. ( |2^ ) is stable, which implies the stability of a trav- 
elling helix solution. However, we can expect that for 
other sets of parameters of the three-dimensional CGLE 
the parameters of Eq. (|2j) may fall into the unstable re- 
gion, e.g the range of amplitude turbulence. This would 
imply chaotic oscillations of the helix. 

We expect that these weakly-nonlinear equations ( [2^ ) 
may also serve as building blocks for the understanding of 
weak vortex turbulence in the CGLE in a certain region 
of parameters. 

An interesting question in this context is: could the he- 
lices form a bound state similarly to spiral waves? One 
could imagine a stable double-helix state, similar to the 
DNA molecule. We have performed a preliminary numer- 
ical investigation of the double-helix configuration. The 
results presented in Fig. ^ indicate that the double helix 
is unstable: the outer helix expands, whereas the inner 
helix shrinks. However, we cannot exlude the possibility 
for existence of a stable double helix in some (narrow) 
parameter range of the CGLE. 




FIG. 11. Sequence of snapshots demonstrating the break- 
down of a double helix for c = 0.5 and e = 0.25. (a) t = 60, 
(b) t = 120 and (c) t = 240. 



VI. CONCLUSION 

We have derived an equation of motion for the vor- 
tex line valid in the high-dispersion limit of the CGLE. 
Using a general linear stability analysis, we have found 
that in a wide range of parameters of the CGLE the vor- 
tex line is unstable with respect to spontaneous stretch- 
ing and bending, resulting in the formation of persistent. 
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dynamic entangled vortex configurations. In fact vortex 
lines in three dimensions are unstable in a much wider 
range of parameters then two-dimensional spiral waves. 
This emphasizes the deficiency of previous approaches 
relating local filament velocity only to local curvature. 

Slightly beyond the onset of the three-dimensional in- 
stability we found stable traveling helix solutions, which 
bifurcate supercritically from the straight vortex. Qual- 
itatively similar helicies have recently been observed ex- 
perimentally in heart tissues [Q. Note that the acceler- 
ation instability in two dimensions is subcritical Q . Un- 
stable helix solutions can maybe serve as building blocks 
for a weakly-nonlinear theory of three-dimensional vortex 
turbulence. 

Let us now discuss implication of our result for the 
well-known phase turbulence problem in CGLE. As was 
found in Ref. ||3l|] , in two dimensions phase turbulence 
is never a global attractor since it is unstable with re- 
spect invasion by defect turbulence. However, in three 
dimensions we may expect that at least in the region 
of parameters away from the stretching instability, the 
vortex rings collapse and the phase turbulence regains 
its stability. But inside the three-dimensional instability 
region we may speculate that there is always a possibil- 
ity for creation of a large enough vortex loop which will 
expand and invade the phase turbulence. 

Our result could be verified in experiments with au- 
tocatalytic chemical reactions in gels in the regime of 
oscillatory instability. The limit of a large dispersion 
b > be can probably be achieved by doping with addi- 
tional chemicals, thus changing the relative mobility of 
reacting components. 

Persistent entangled vortex configurations are known 
from numerical simulations of excitable reaction-diffusion 
systems p^JT^ . Our preliminary investigation of 
reaction-diffusion systems shows that the under-damped 
core dynamics here is also responsible for long-lived vor- 
tex loops and persistent entangled vortex configuration 
p2| . In this case, the expansion of the vortex loops is 
not necessarily related to a " negative line tension" of the 
filament, but again is the manifestation of acceleration 
effects, similar to the situation in the CGLE. 

Recently, the amplitude equation governing the dy- 
namics of an elastic rod was derived p3| . The structure 
of solutions, found in Ref. [ ^3t are remarkably similar to 
those of the CGLE. It is plausible to assume that, in 
some distinct range of the parameters, the equations of 
motion of the twisted elastic rode can be reduced to the 
equations for the vortex line in CGLE. We also speculate 
that our results are relevant for inviscid hydrodynamics. 
In the limit of fe, c — > cxd, Eq. (|l|) reduces to the defo- 
cusing nonlinear Schrodinger equation (NSE), which is a 
paradigm model for compressible inviscid hydrodynam- 
ics. Although the vortex lines are stable in the framework 
of the NSE, the corrections arising from the CGLE cause 
their destabilisation and stretching. 
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APPENDIX A: LIMIT OF NONLINEAR 
SCHRODINGER EQUATION 

1. Adjoint mode 

In order to derive the friction coefficient x in Eq. ( p^ ) 
we have to fulfill the solvability condition. In Sec. II this 
was done numerically. The solvability condition means 
the orthogonality of the r.h.s. of Eq. ( JTs] ) to the adjoint 
mode of Eq. (p5|). The adjoint equations are of the form: 



-2F^B+ = 



1 ~^ \ Qj. 



'-^{r^'FAt)-'-^ 



= 0. 



(Al) 



The solvability condition of Eqs. (|l^) can be expressed 
in terms of functions A'^ , Bf : 



dtvh - 2evl2 = 0, 



(A2) 



where 
h 



r^drFAt 



h = 



rdr {F'A+ + [F^' + iF/r)B+) 



(A3) 



From Eq. (A2) we readily obtain the friction coefficient 
(compare with Eq. (ph) 



X 



-2/2//1 



(A4) 



For arbitrary c the localized adjoint mode A^, B^ can 
be determined only numerically 1^. However, for c — > 
00, i.e. in the limit of the perturbed NSE, the adjoint 
mode, and can be calculated fully analytically. For 
c — > cx) the selected wavenumber fco vanishes (see, e.g. 
p5[), the solution approaches the vortex solution of NSE, 
and Eq. (Al) becomes self-adjoint. We introduce a small 
para meter — 1/c ^ 1 and assume /i 3> e. Equations 
(^) read: 
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/ F2 

AA+ ~ 2 — A 

V M 

+ 2F^B+ = 



(A5) 



For ji — Eqs. ( [A5| ) are self-adjoint and the localized ad- 
joint eigenmode conincides with the complex-conjugated 
translation mode 



At 
Bt 



F' 

-iF/r 



(A6) 



In order to evaluate Eqs. (A4) we take into account that 
for any /i 7^ the functions , B ^' decay exponentially 
for r —^ 00. Therefore Eq. (A6) can be considered an 
approximation valid within a finite interval < r < Rq, 
where the cutoff i?o 3> 1 will be determined from the 
matching condition with the outer asymptotics of the 
solution. For r — > cx3 we can simplify Eqs. (A5) using 
that ^ 1 - ^/r^ + and fco < 1: 

kAt - 2 (1^+ + l^ir^'Bt) 2Bt = 



^Bt+2{l§-^iri.'At)-'-^]=0. 



(A7) 



From the first Eq. (A7) we can explicitly express A^ in 
terms of B^ , because for r — + cx3 all terms in the first Eq. 
(A7) except A^ / IJ. and B^ ■ In the first relevant order we 
obtain 



At = ^^Bt 



(A8) 



Substituting now Eq. (A8) into the second Eq. (A7), 
and dropping higher-order terms one has 



d^M + -drB+ \b+ + ^^ir^P'Bt) 
r r'' r or 



0. (A9) 



Eq. (A9) is reduced to Bessel's equation by the substi- 
tution Bt = S cxp{— flip) , leading to 

d^S +^drS - in^k^ + = 0. (AlO) 

It has the localized solution S ~ Ki{fikor). Thus, the 
outer solution is of the form: 



At 
Bt 



C{^^]eM-y^)Ki{fikor) (All) 



where C is a constant determined from the matching with 
the inner solution Eq. (A6). For r — > cx) the functions de- 



cay exponentially: At , B ^ exp[— 2/i|fco |r]/r^/^, smce 
ijj kor. The solutions (A6) and ( All ) match in the 



intermediate region r 3> 1 and /ifepr <C 1. Expanding 
Ki{fikQr) for small arguments, we find C = —ifikQ. 



2. Friction coefficient 



To calculate the integrals (A3) we introduce the cutoff 
i?o = l/(A'fco) ^ 1 and split the interval of integration 
into two parts: < r < Rq and Rq < r < 00. In 
the first interval we us e the inner representation of the 
eigenf unctions Eq. ( |A6|), and in the second one the outer 
representation Eq. ( |A11 ). For the inner interval we ob- 
tain 



11^ 



= / rdr {{F'f + {F/rf) =h+ logi?o, (A12) 
Jo 

where Bq and 61 are some constants. From the outer 
integration for I2 w e have (using —> k,F —> \ for 
r ^ 00, see Eq, (A3) ) 



r^drFF' = 60 +A*logi?o 



1° w / r{ko+i/r)drBt 



Ro 
0.884 



log(i?oAi^o) ~ j0.666//i. 



(A13) 



In order to evaluate the integral Ii we need the next order 
of the function At to compensate the log arit hmic diver- 
gence of the inner integral. From Eq. (A5) we obtain 
At = fi{Bt + iBt/r^)- Thus, we have 



r^dr[—i^ko{l + — )]/iii'i(/ifcor) exp(^A:or) 



-0.884/x- /xlog(i?oAi^o) - — 



OAi 
fik^ 



(A14) 



Combining now the outer and inner expansions, we see 
that Rq drops out. The friction coefficient is of the form: 



X 



J - log(^fco) + ci - iO.666/^ 
'-^log(^/co) + Co - iOA/{fik^) ' 



(A15) 



where ci = 61 - 0.884, cq = fco - 0.884^. Now, us- 
ing Hagan's expression for the selected wavenumber p^ ] 
in our scaling of the CGLE parameters one finds kQ w 



2/1" 



-3/2 



exp[-7r/(2/i) 
13.3 



■ exp 



7" 

TT 



3.45 



■ exp 



0.098]. We finally obtain 
- 27-0.196^ 

(l + 3.92i). (A16) 



1 



ITT 

08 



The real part of the friction coefficient is always nega- 
tive, which implies the instability of the spiral core in 
two dimensions and ia stretching instability in three di- 



mensions. 
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APPENDIX B: RELATION BETWEEN 
ACCELERATION INSTABILITY AND 
STABILITY PROBLEM 

The result of Sec. II for weakly-curved vortex filaments 
can be formally derived from the linear stability problem 
Eq. through systematic expansion in e up to second 
order. The eigenvalue problem for e <^ fi is oi the form: 



[1] 
[2] 
[3] 
[4] 
[5] 



) ^ (Ae + + XBi [6] 



fi or t 

or r 



Now we expand Eq. (Bl) in e. The solution is repre- 
sented in the form 



Bl) ^ \Bi^) 
A = eA(i) - 



^(1) 



^a(2). 



(B2) 



We consider k ^ 0(e) and denote k — kjt. 

At zeroth order in e we simply obtain B^'')) = 



0, where L is the r.h.s. ofEqs. (Bl). Clearly the solution 
is the translation mode. 



F' 



(B3) 



At first order in e we obtain 
Li 



Bwj~ y -F' J- 



The equations have an exact solution corresonding to the 
'family mode', which exists for e = 0: 



AW 



aw / 

2 [rF 



At second order in e we obtain equations 
^(2)^ /(XW +k2)F' - ^^rF^ 



L 



B( 



(Ad , . 
(A(i)+fc2)(F^' + i£) 



A(2), 



(B5) 



(B6) 



These equations have a bounded solution if the solvability 
condition is satisfied. It is easy to see that Eqs. ( p36D are 
identical to Eqs. ( |T^ ) if we take into account w ~ A, dtv — 
Xv and k — k^. The solvability condition implies that 



(A(i))2 + x(A(i) + P) = 0. The last term in Eqs. (BE) 



can be omitted, since it gener ates a nonsingular solution 
(family mode, compare Eq. (B5)). Thus, we reproduce 
the result for a weakly-curved vortex Eq. (|l^) . 
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